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Abstract 

'Umbral calculus' deals with representations of the canonical commutation relations. 
We present a short exposition of it and discuss how this calculus can be used to 
discretize continuum models and to construct representations of Lie algebras on a 
lattice. Related ideas appeared in recent publications and we show that the examples 
treated there are special cases of umbral calculus. This observation then suggests 
various generalizations of these examples. A special umbral representation of the 
canonical commutation relations given in terms of the position and momentum 
operator on a lattice is investigated in detail. 



1 Introduction 

Umbral calculu^ is an analysis of certain representations of the commutation relations 

[Qi , Xj] = Sij 1 , [Qi , Qj] = = [xj , Xj] (1.1) 

in terms of operators on the algebra of polynomials in variables Xi, i = 1, . . . ,n (see 0] 
for reviews). In particular, it provides us with representations by operators acting on 
polynomials of discrete variables. Let us assume that Qi, Xj is such a representation and 
let A{yi, d/dyj) f{yk) = be a differential equation on IR" with a polynomial solution f.f\ 
Introducing multiplication operators y^, we can write it in the form ^ 

A{yi,d/dy,) /(yfc) 1 = 0. (1.2) 

^This terminology goes back to the nineteenth century mathematician Sylvester who used the Latin 
word umbra to denote something which would nowadays be called a linear functional. See also ^ . 
^Here and in the following an expression like f{yk) stands for /(yi, . . . , y„). 

■^In this expression the 1 plays the role of a 'state' on which we act with an operator algebra to generate 
an irreducible representation of the latter. 



The operators and d/dyj do satisfy the commutation relations ( |l.lj ), of course. The 
verification that f{yk) solves the original differential equation is now translated into an 
algebraic problem which only requires the abstract commutation relations i.e., it 

does not depend on the specific choice of representation. Defining f{xk) '■= fi^k) 1, then 
also 

A(x,,Q,)/(x,) =0 (1.3) 
holds which is a difference equation. We have simply substituted 

d 

^ Xi , — ^ Qj . (1.4) 

If f{yk) solves the original differential equation, then f{xk) is a solution of the correspond- 
ing difference equation. 

For differential equations possessing polynomial solutions, the notion of quasi-exact 
solvability has been introduced ||^. Several examples are provided by eigenvalue problems 
in quantum mechanics. A corresponding example for the above discretization procedure 
appeared recently in 0. In section ^ we show that its underlying structure is umbral 
calculus. 

The above operator substitution yields a mapping of an eigenvalue equation for a 
differential operator to an eigenvalue equation for a difference operator together with a 
'formal' mapping of solutions. It seems that we have a general procedure for 'isospectral 
discretization' of differential operator eigenvalue problems. The problem, however, is that 
(besides for polynomials) the mapping of solutions in general only works at the level of 
formal power series, but does not respect convergence properties. 

Also, in the abovementioned treatment |Q of eigenvalue problems one does not re- 
ally get a discretization of the original quantum mechanical problem since that involves 
non-polynomial functions. For serious applications we therefore need an extension of the 
procedure sketched above beyond polynomials and formal power series. Such a discretiza- 
tion method could then be of interest for solving differential equations numerically. 

The commutation relations of differential operators A(yj, d/dyj) and -B(yi, d/dyj) are 
preserved under the substitution ( |1.4| ). In this way one obtains representations of operator 
algebras, in particular Lie and Hopf algebras, by operators acting on functions on a lattice. 
An example appeared recently in where representations of the Poincare and the k- 
deformed Poincare algebra on a lattice were constructed. In section |] we explain how 
it fits into the umbral framework. 

All this raises the question whether it is possible to understand (some of the) umbral 



maps ( p..4|) on algebras of non-polynomial functions. In view of possible applications to 
quantum mechanics, it would be of interest to have Qi,^j defined on the Hilbert space 
of square summable functions on a lattice. Is it possible that Xj and —iQj (which as a 
consequence of ( |1 . 1| ) satisfy the canonical commutation relations of quantum mechanics) 
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are selfadjoint operators and is ( |1.4|) perhaps a unitary equivalence ? Our work intends 
to contribute to the clarification of such questions. An example of particular interest is 
suggested by the work in 0. The representation of the canonical commutation relations 
which appeared there is investigated in detail in section 

Section ^ contains a brief introduction to our understanding of umbral calculus. By 
no means it intends to cover the whole subject. An example treating symmetries on a 
lattice is then presented in section |[ Another application is discussed in section H, partly 
motivated by In section ^ we slightly generalize the umbral framework of section |^. 
We also comment on a representation of the Poincare algebra on a lattice which appeared 
in 0]. Its underlying representation of the canonical commutation relations is the subject 
of section ^ It leads us to a framework for quantum mechanics on a lattice. Some 
conclusions are collected in section 0. 

2 A brief introduction to umbral calculus 

In this section we recall some notions and results from umbral calculus. We refer to 0, 
for the corresponding proofs and further results. For simplicity, we restrict our consider- 
ations to the case of a single 'coordinate' x. All results extend to several (commuting) 
variables in an obvious way. 

An operator O acting on the algebra (over a field of characteristic zero, like R or C) 
of polynomials in x is shift-invariant if it commutes (for all a in the field) with the shift 
operators Sa (defined by Saf{x) = f{x + a)). 

The Pincherle derivative of an operator O is defined as the commutator 

O' ■.= [0,^]=0x-^0 (2.1) 

where x is the multiplication operator, acting on polynomials in x by multiplication with 
X. The Pincherle derivative of a shift-invariant operator is again a shift-invariant operator. 
The umbral algebra is the algebra of all shift-invariant operators. The Pincherle derivative 
is a derivation of the umbral algebra. 

A delta operator Q is a linear operator, acting on the algebra of polynomials in x, 
which is shift-invariant and for which Qx is a nonzero constant. It can be shown that 
Q'~^ exists (as a linear operator on the space of polynomials) and commutes with Q. If 
we define 

x:=xg'~^ (2.2) 

it follows that 

[g , x] = I (2.3) 

where I stands for the identity operator. In this way each delta operator Q provides us 
with a representation of the canonical commutation relations on the algebra of polynomials 
in X. 
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A polynomial sequence qk{x), k = 0,1,2 is a sequence of polynomials where qk{x) 
is of degree k. A polynomial sequence is called basic for a delta operator Q if go(^) = !> 
qk{0) = whenever k > 0, and 

Qqk = kqk^i . (2.4) 
It turns out that basic sequences are of binomial type, i.e., they satisfy 

qk{x + y) = Yl (^^ *(^) 1k-e{y) ■ (2-5) 

The basic polynomial sequence for Q is given by 

gfc(x) =xgfe_i(x) =xM (2.6) 

which is known as the Rodrigues formula. 

An operator which maps a basic polynomial sequence into another basic polynomial 
sequence is called an umbral operator (0], p. 28). Defining 

/» := /(x) 1 (2.7) 

for a polynomial /, ( |2.6|) shows that the operator ~ is an umbral operator. 
An associative and commutative product is defined by 

/>)*/i(x):=/(x)/i(x)l. (2.8) 

In particular, qk{x) * qi{x) = qk+i{x). The delta operator Q is a derivation with respect 
to the *-product, i.e., 

Q[p{x) * q{x)] = {Qp{x)) * q{x) + p{x) * Qq{x) (2.9) 

for polynomials p and q. 

Example 1. For Q = d/dx we have Q' = 1 and therefore qk{x) = x^ which is the simplest 
polynomial sequence. 

Example 2. Let Q = D/{D -1) with D := d/dx. Then Q' = -{D - l)-^ and qk{x) = 
[—X (D — 1)^]^ 1 are the basic Laguerre polynomials P]. D 

As pointed out in the introduction, we are particularly interested in the case where the 
algebra of polynomials in x can be realized as an algebra of functions on a discrete set. In 
the following two examples we may choose x to be the canonical coordinate function on 
an infinite lattice with spacings a (where a is a positive real number). In the way outlined 
in the introduction, both examples provide us with a prescription to translate functions 
on IR and differential operators into corresponding functions and operators on a lattice. 
The interesting aspect is that this prescription not only maps a differential equation into 
a corresponding difference equation, but it also allows us, in principle, to calculate the 
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solutions of the difference equation from those of the differential equation (see sections ^ 
and H). 

Example 3. Let Q = where d+ is the forward discrete derivative operator, 

(9+/)(x) = -[/(x + a)-/(x)] (2.10) 
a 

acting on a function /. We find {Q'f){x) = f{x + a) and therefore Q' = Sa, the shift 
operator. Hence, 



(x^-i)^ l = x{x-a)---{x-{k-l)a) = x^^^ (2.11) 



where x*^'^-' is the fc-th (falling) factorial function^ Some formulas for the ^-product associ- 
ated with the discrete derivative delta operator can be found in the appendix. Analogous 
formulas hold for the backward discrete derivative operator (9_ which is formally obtained 
from (9+ replacing a by —a. D 

Example 4- For the central difference operator |^ 

Q = ^iS,-S-a) = l{d+ + d^) (2.12) 

we have Qf{x) = [f{x + a) — f{x — a)]/(2a). Solving (|2.4| ), one finds the basic sequence 

fc-i 

qk{x) = X Y[{x + ka - 2na) {k > 1) (2.13) 

n=l 

and qo{x) = 1, qi{x) = x. Furthermore, 

Q' = ^iSa + S.,), Q" = a^Q. (2.14) 
Using the Rodrigues formula, 

k 

Q'-^qk{x) = \l\x + {k + l)a-2na] (fc > 0) (2.15) 

n=l 

which shows that Q'~^ is indeed well-defined on polynomials in x. The operator Q'~^ also 
exists as a selfadjoint operator in the Hilbert space £2(0-^), see section ||. D 

Example 5. Over a finite field there are finite-dimensional representations of the commu- 
tation relation (p^.3|). Over the matrices 
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^For a — 1 and x S IN it counts the number of injective maps from a set of n elements to a set of x 
elements. 
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provide us with an example which generahzes in an obvious way to the Galois fields GF{p"') 
(where p is a prime and n G IN). Though in this case we leave the usual umbral framework 
since we consider a field which is not of characteristic zero, some basic constructions and 
results remain valid. D 

As long as we restrict our considerations to operators acting on polynomials, everything 
works smoothly. We are, however, also interested in more general classes of functions 
and in particular power series. In general, an umbral operator like ~ does not preserve 
convergence of such a series. The result of the application of an umbral operator to a 
power series a priori only makes sense as a formal power series. A problem is then to 
determine its domain of convergence (which may be empty) and a possible continuation. 
It seems that little is known about the convergence of power series obtained via umbral 
maps. 



3 Symmetry operators on a lattice: an example 

In this section we generalize the example 3 of section ^ to n dimensions. As an application 
of the umbral method, a representation of the Lie algebra of 5*0 (3) on a lattice is then 
presented. Let be the canonical coordinate functions on an n-dimensional 

(hypercubic) lattice with spacings a^. We define delta operators 

iQif){x) ■= — [f{xi, Xi-i, Xi + a,, Xi+i, . . . , x„) - f{x)] (3.1) 

acting on functions of ). The corresponding Pincherle derivatives are the 

shift operators Si acting on functions as follows, 

{Sif){x) := /(xi, . . . , Xi_i, Xi + ai, x^+i, . . . , a;„) . (3.2) 

The operators Xj = Xj S^^ and Qj then satisfy the commutation relations ( |1 . 1| ) on the 
algebra of polynomials in the variables xi, . . . , Xn- 

As outlined in the introduction, given a representation of a Lie algebra in terms of the 
operators and d/dyj acting on functions on IR", ( |1.4D maps it into a representation by 
operators acting on functions on a lattice. For the angular momentum operators in three 
dimensions this means 



d 

Li = -'^Yl ^ijk Yj TT- ^ Li := -i Eijk Xj Qk (3.3) 

j,k j,k 



where 



Lif{x) = -i '^eijkXj {Qkf){x - aj) (3.4) 



using the notation x — aj = {xi, . . . , xj-i, xj — aj, Xj+i, . . . , a;„). 
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What are the corresponding 'spherically symmetric' functions on the lattice ? We have 
to find the solutions of Lj/(x) = 0. From the corresponding solution in the continuum 
case, we know that / should depend on only through Yl\=i x^l = Y%=i Xk{xk — a^) 
or ^-products of this expression. The set of lattice points determined by the equation 
J2k=i ^ki^k — cik) = constant therefore constitutes the analogue of the 2-sphere in the 
continuum case. Of course, only for special values of the constant it will be non-empty. 
For a lattice with equal spacings in all dimensions, the mappings Xk ^ xe and Xk ^ a — Xk 
leave the above expression invariant and thus help to construct the 'lattice spheres'. 

4 Isospectral discretization of eigenvalue equations 
via umbral calculus ? 

In differential equations were called 'quasi-exactly solvable' if there is at least one 
polynomial solution and 'exactly solvable' if there is a complete set of polynomial solutions. 
The relevance for physics has been established in a series of papers where quantum 
mechanical eigenvalue problems were collected which can be reduced to equations having 
polynomial eigenf unctions via an ansatz of the form 

m = mf{y) (4.1) 

with a fixed non-polynomial function / on IR. The most familiar example is provided by 
the (one-dimensional) harmonic oscillator. In this case 

^{y) = <P{y)e~y'" (4.2) 

converts the Schrodinger equation into a differential equation for which has the Hermite 
polynomials as a complete set of solutions. Another example is the radial part of the 
Schrodinger equation for a hydrogen atom. 

In a discretization procedure has been proposed for a differential operator eigenvalue 
equation possessing polynomial solutions such that the resulting difference equation has 
the same spectrum. It corresponds to an umbral map in the sense of section |^ with 
the choice Q = <9+, the forward discrete derivative operator.]^ However, the procedure 
does not work well, in general, when applied to the original eigenvalue problem which we 
started with. Though we do get a discrete eigenvalue problem in this way which is naively^ 
isospectral, problems arise when we try to translate the non-polynomial solutions. This 
will be illustrated with the following examples. 

^The umbral framework provides us with several alternatives, of course, which have not been considered 
ini. 

^In order to formulate a well-defined eigenvalue problem, we have to specify a suitable function space 
in which we are looking for solutions. Each eigenfunction of a differential operator is mapped to an 
eigenfunction of the corresponding discrete operator (or at least a formal power series which satisfies the 
discrete eigenvalue equation). Note, however, that the discrete equation may have additional solutions. 
In particular, this is the origin of boson or fermion doubling in lattice field theories (cf [|[ 0). 



7 



Example 1. We apply the umbral map to a simple differential equation, 

^f{y) = kf{y) ^ Qf{x) = kf{x). (4.3) 
dy 

From the solution f{y) = exp{ky) of the differential equation the corresponding solution 
of the difference equation on the rhs of ([4 .31) is then obtained as follows, 



/(x) = /(x)l = ;^-xW. (4.4) 
e=o ^■ 

Though we would like to choose x as the canonical coordinate function on the lattice aZ, 
it may be helpful at this point to consider it as a coordinate function on IR in view of a 
possible analytic continuation of the power series obtained from the umbral procedure. A 
priori, we obtain / only as a formal power series. For real k, the series in ( [4.4|) (which is 
a special case of a Newton series) converges everywhere on the real line if \ka\ < 1. For 
\ka\ > 1 the series is everywhere divergent, except for non-negative integer multiples of a 
(see |TT]], for example). The difference equation on the rhs of ( |4.3| ) has no nonvanishing 
solution for k = —1/a. For all other values of /c G C the solutions are given by 

/(na) = /(0)(l + A;a)". (4.5) 

With /(O) = 1 this extends the series obtained above (for k > —1/a). D 

Example 2. Let us now apply the umbral map to the Hamiltonian of the harmonic 
oscillator, 

H = --^ + -y^ ^ H = --Q'' + -kK (4.6) 

2 dy^ 2^ 2^ 2 ^ ^ 

The eigenvalue equation for H is then translated into the following eigenvalue equation 
for H, 

Hi){x) = ^ [-Q^{x) +x{x - a)ij{x-2a)] = E i}{x) (4.7) 

which is a difference equation {Q = d+). From the solution ipo{y) = exp(— y^/2) of the 
original eigenvalue problem we obtain the solution 

1=0 ^ ^■ 

of the difference equation ([4.7| ) with = 1/2 as a formal power series. Using 

a;(2^+2) _{x-2ia){x-2ia-a) x^^^^ 
2^+1 (£+1)! ~ 2(£ + l) 2n\ ^ '^^ 
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the quotient criterium shows that the series is everywhere divergent, except at values of 
X which are non- negative integer multiples of a (where the series terminates). ([4 .81 ) thus 
only determines a solution of the difference equation 

x{x-a) tjj{x - 2a) = \[ip{x) -2i){x + a)+ ipi^x + 2a)] + ijj{x) (4.10) 

on the non-negative part of aZ. The Ihs of (|4.10|) vanishes for x = and x = a. Also 
the rhs vanishes if we calculate the corresponding values of ip using ([4.8|) . Our solution 
can therefore be extended to the whole of aZ. But the extension is not unique since 
ipl—a) and ip{—2a) can be chosen arbitrarily. This shows that the difference equation has 
more independent solutions than the differential equation we started with. The umbral- 
mapping of ipo can, however, be completed to yield a solution of the difference equation 
which exists everywhere on aZ. This is done by expanding i/jq into power series about 
negative multiples of a and acting with ~ on these series. 

The higher eigenfunctions of the harmonic oscillator are products of Hermite polyno- 
mials with ipQ, 

My)='Hn{y)My) ■ (4.ii) 

Now ipn{x) is obtained by replacing the ordinary product by the *-product (cf the ap- 
pendix), ipo{y) by ipo{x) as given above, and the Hermite polynomials by the 'discrete 
Hermite polynomials'. The latter are obtained from the generating function 

e=o ^- e=Ok=l '^r 

hy Unix) = {d/ds)^F{x,s)\s=^. □ 



5 Some more umbral calculus 

There is a generalization of the calculus described in section |^. Given a representation 
of the commutation relation ( |2.3| ) by operators Q and x as in section |^, and given an 
operator A on the space of polynomials which commutes with Q, then the new operator 
X + A together with Q also satisfies the commutation relation. In the following, let x 
denote such a more general choice (than the special one in (|2.2|) ). Defining 

Sk{x):=yiH (5.1) 

one finds 

Qsk = ksk-i (5.2) 

by use of the commutation relation (|2.3| ). Such a polynomial sequence Sk is called a Sheffer 
set for the delta operator Q in the umbral literature. The basic polynomial sequence qk 
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for a delta operator Q is a special Slieffer set. If Sk is a Sheffer set for Q, then there is 
an invertible shift-invariant operator which maps the Sheffer polynomials to the basic 
polynomials g^. Furthermore, 



s„{x) = J2 ylj'^k{0)qn-kix) . (5.3) 

We refer to for proofs and further results. Again, we define f{x) := /(x) 1. 

We have to stress that not all umbral results established for the special choice (|2.2|) 
for X translate to the more general case considered in this section. In general, Sfc(O) ^ 
and the Sk are not binomial. 

A particularly interesting choice for x turns out to be 

x = i(xg'-^ + Q'-^x). (5.4) 

From umbral calculus we know that Q'^^ commutes with Q. One can then easily verify 
that [Q,x] = I. The advantage of ( ^.4|) over (|2.2|) is that it is more symmetric and thus 
opens the chance to turn x and iQ into Hermitian operators on a Hilbert space.[] For 
Q = d/dx we have Sk{x) = qk{x). In case of the (forward) discrete derivative operator 
one finds 

1 13 9k — 1 

Sk{x) = ^(x5-^ + S^'^r l = {x--a){x-^a)---{x- ^— a) . (5.5) 



Another realization of ( p. 41 ), involving the central difference operator, will be the subject 
of the following examples and the next section. In that case, we have 



So{x) = l, si{x)=x, S2{x)=x^- — , ... (5.6) 
using i JU^ ) and Ql = 0. 

Example 1. Let us consider again the example of the harmonic oscillator. Using ( ^.4|) and 
the central difference operator ( p.l2| ), the corresponding Schrodinger equation is umbral- 
mapped to 

i |v5(x) = \ + g'-^(x^ - ^) + 2a' Q'-'Q^ + ^-a' Q'-'Q'] ^{x) (5.7) 

where on the rhs we have naively commuted all the non-local operators Q'^^ to the left. 
Acting with Q"^ on this equation results in a finite difference equation (with respect to the 
space coordinates). However, if we discretize the timeQ in order to solve the initial value 

X is Hermitian and Q anti-Hermitian, then Q' and Q' ^ are Hermitian and thus also the operator 
in (O ). 

**This can be achieved via an umbral map, of course. 
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problem for the above equation on a computer, calculation of the wave function at the 
next time step requires Q'~^. But to explore an equation of the type above numerically, we 
have to use an approximation with a finite lattice. Choosing periodic boundary conditions 
(i.e., a periodic lattice), there are convenient formulas for Q'~^. On a periodic lattice with 
N = 2m sites where m is odd, the equation 

m— 1 

Q'-'=Y^{-lfSf^\ (5.8) 

fc=0 

holds.0 For odd N one finds instead 

(N-l)/2 {N-l)/2-l 

Q'-'= J2 (-l)^+(^-i)/2 5f + (-1)'^''+'- (5.9) 

A;=0 fc=0 

In the following section, the quantum mechanical setting behind ( [5.7| ) is investigated more 
rigorously. D 

Example 2. In section |] we determined the 'lattice spheres' with respect to some umbral 
representation. Instead of ( |2.2| ) here we choose 

X, := X, {Si + Sry^ + {S, + Sry^ x, (5.10) 

with the shift operators defined in (|3.2| ). This means that we consider (|5.4]) generalized 
to several dimensions with central difference operators 

Q. = ^(^.-^r')- (5.11) 

Using Q'- = afQi and Qil = 0, we find the following equations for 'lattice spheres' in 
three dimensions, 

f](x,)2 1 = - f ] = constant . (5.12) 

A spherically symmetric potential on the lattice is then a function which only depends on 
^-products of Y.\=A^Xif - a?/2]. D 

Example 3. A familiar representation of the Poincare algebra is 

d 

-P/x = -i^, Mi = 'Yeijkyj Pk , Li = yoPi - K-yiPo . (5.13) 
•^^M j,k 

These operators act on functions on IR^ (with canonical coordinates y^). The commutation 
relations are then preserved when we perform in the expressions ( |5.13|) the substitutions 

d 

ji ^ Xi , — Qi , 1 = 1,2,3 (5.14) 

oyi 



^For even m, Q' is not invertible. 
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with the operators defined in (^.10|) and ( p.llj ). In this way we obtain a representation 



of the Poincare algebra on a lattice with spacings Oj (and continuous time as long as 
and Pq £^re kept unchanged) .P^ The quadratic Casimir operator of the Poincare algebra 
in this representation is 

c = -d',+Y.Ql- (5.15) 

k 

where dt := d/dxQ. There is, however, a drawback of the representation presented above 
and also those given in 0. As pointed out in 0, the Klein-Gordon equation built with 
the operator (|5.15|) suffers from a boson doubling problem analogous to the more familiar 



fermion doubling problem in lattice field theory (see |10|, for example). This leaves us 
with a Poincare-invariant theory with 8 species of bosons. If the time dimension is also 
discretized, one obtains 16 species. The relation between this Klein-Gordon equation and 
the Dirac equation for 'naive lattice fermions' is the same as in the continuum, 

(i 7°5t + i - m) (i -f'dt + i -f'Qk + m) = -d^ + ^ _ ^2 _ ^^_^q^ 

k 

The representation of the Poincare algebra acting on continuum spinor fields is mapped via 
( ^.14| ) to a representation on the lattice which leaves the lattice Dirac equation invariant. 

□ 



6 Via umbral calculus to quantum mechanics on a 
lattice 

In this section we investigate the umbral discretization method with the central difference 
operator Q = {Sa — S^a) / (2a) and the symmetric operator ( ^.4|) . It will be shown that they 
define selfadjoint operators on the Hilbert space £2(0..^), the space of square summable 
functions on the infinite lattice with spacings a. We thus have a rigorous framework to 
explore the 'umbral map'. 

By standard arguments x is selfadjoint with domain {/ G i2{aZ) \ x/ G £2{aZ)}. 
The Fourier transformation f v-^ F where 

f[x) = ^ r''' F{k)(^^''dk (6.1) 

V27r J-ir/a 



is an isomorphism £2(a^) — > L'^ia- Here and in the following Li stands for L^([— 6, 6]), 
the space of square-integrable functions on the interval [—6, h] . It is more convenient for 
us to define the domain of x now as follows, 

TT TT dF 

V^ = {f e i2(a^) I F absol. continuous, F( ) = F(-), — G LLJ . (6.2) 

dk ' 



a a 
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See also [||. A representation on a four-dimensional space-time lattice is obtained by extending the 



map (5.14) to yo and d/dyo. 
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The action of x on £2(0-^) then corresponds to the action of id/dk on the domain in L^/^ 
specified above.[^ Its spectrum is {raa | G Z}. 

Next we note that —iQ is a bounded selfadjoint operator on £2(0.^)- In -^^/^ it acts by 
multiplication with sm{ak)/a. Concerning the umbral map we can conclude the following, 

operator: —id/dy 1— >■ — iQ 

spectrum: R { A G R | |A| < i} 

eigenfunctions: fxijj) = exp(iA?/) f\{x) = exp[i(a;/a)arcsin(Aa)] 

The eigenfunctions of —iQ can indeed be calculated directly from the power series ex- 
pansions for those of —id/dy (with the help of 0, section 6.5). The spectrum of — iQ is 
bounded, however, in accordance with the boundedness of the operator. Only in the limit 
a — s> we recover the full spectrum of the continuum momentum operator. 

For / G Px we have Qf G Px- The operator Q' = [Q, x] = (Sa + S-a)/'^ is then defined 
on Dx- It is bounded and can be extended to a selfadjoint operator on i2{aZ). Q'f = 
for / G £2{aZ) implies / = 0. Hence Q'~^ exists on "Pq'-i = Q'(£2(a-^)) and is selfadjoint 
(Lemma XII. 1.6 in |T^). The Fourier transform of Q' acts in L"^/^ by multiplication with 
cos(aA;). The operator Q'~^ therefore acts by multiplication with 1/ cos(aA;) on the domain 
{FeLl/jF{k)/cos{ak)eLl/J. 

It remains to investigate the operator x = (xQ'^^ + Q'^^x)/2 which is Hermitian on 
the dense domain 



v^ = {fev^n Vq,-. I x/ G Vq,-., Q'-'f g v^} . (6.3) 

Without any calculations we can immediately conclude the following. The operator x^ can 
be defined on a dense domain on which it commutes with complex conjugation. Accord- 



ing to Theorem XII. 4. 18 and Corollary XII. 4. 13(a) in [12| this operator has selfadjoint 



extensions. Let us recall a theorem due to Rellich and Dixmier (see Theorem 4.6.1 in 



0) 



Theorem. Let q and p be closed Hermitian operators on a Hilbert space Ti. such that 

(1) [p, q] = — i on a subset Q C Pq fl Vp dense in Ti which is invariant under q and p, 

(2) p^ + q^ on n is essentially selfadjoint. 

Then p and q are selfadjoint and unitarily equivalent to a direct sum of Schrodinger 
representations. 

An isomorphism i2{aZ) = L^(R) maps the operators x and — iQ to corresponding opera- 
tors in L^(R). These operators cannot be unitarily equivalent to those of the Schrodinger 
representation since the latter are both unbounded. Besides (2), the operators x (which 



has a closed Hermitian extension |T2[) and — iQ fulfil all assumptions of the last theorem. 

^^The latter is a standard textbook example of a selfadjoint operator. Via Fourier transformation it is 
mapped to a selfadjoint operator on C l2{a,'Z\ 



13 



Taking into account that {—iQ)^ is selfadjoint and bounded, it follows that is not essen- 
tially selfadjoint. Together with our previous result this means that has inequivalent 
selfadjoint extensions. 

We now turn to a closer inspection of the operator x which, via Fourier transformation, 
is translated into the operator 

^i/l(i^(il\ (Q 4:) 

2 \cos{ak) dk dkcos{ak)J 

with domain Vx C L'^/a determined by (O). This operator is singular at k = ±7r/(2a) 
and functions in Vx vanish at these points. Assuming that the latter also holds for func- 
tions in the domain of selfadjoint extensions]^, the eigenvalue problem for X separates 
into two independent eigenvalue problems, namely for the following two operators. 

(a) X(i) is X restricted to Vxw := {F G Ll^^^ \ F{k)/cos{ak) abs. cont., F(-7r/(2a)) = 
= Fin/i2a)),XFeLl/,J 

(b) is X restricted to := ^ | F{k)/cos{ak) abs. cont., F(-7r/(2a)) = 
= F(7r/(2a)),F(-7r/a) = F{7r/a),XF e L^} where := L'^ {[-ti / a, -n / {2a)] U 
[7r/(2a),7r/a]). 

In both cases we perform a change of coordinate 

p = - sin(afc) . (6.5) 
a 

Then, with the separation 



F{k) = ^\cosiak)\x{p) (6.6) 

we find for i = 1,2, 



(X Wf) (k) = i ^\cos{ak)\ ^xip) . (6.7) 
The two operators X^^^ now both translate into the more familiar one 



d r2ii /1\ dx 

- on {xGL,/Jxabs. cont., x(--)=0 = x(-), ^ 



Let and F^2) denote the restrictions of F G T>x to [—Tr/a, Tc/a] and [— vr/a, — tt/ (2a)] U 
[tt/ (2a), 7r/a], respectively. Then 



(F, F') = Fiky F'{k) dk = (x(i), X(i)) + (X(2), X'(2)) (6.9) 

J —■K/a 



^^An inspection of the adjoint of X (which contains all selfadjoint extensions) indicates that this holds 
indeed for all selfadjoint extensions. For our purposes it is sufficient to verify a posteriori that it holds 
for all the selfadjoint extensions which we construct below. 
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where 



(XW, X'w) = /'^" xUip) dp {i = l,2). (6.10) 



The selfadjoint extensions of the operator ( |6.8|) are known to be given by 
Da = i^ on := {x G L?/^ | x abs. cont., x(^) = e^^^" x(-^), ^ G (6.11) 
where a G [0, 1) ||l4l. A complete orthonormal set of eigenfunctions of is 



X^(p) := exp[-i(a + n)7rap], nG^, (6.12) 

and Da has a pure point spectrum {{a + n) it a \ n G Z}. A selfadjoint extension of the 
operator x is now obtained by choosing any pair from the set of operators D^- It then 
defines operators X^^^ and Xj^^ and in this way a selfadjoint extension Xq,^ of ^ with 
spectrum {{ai + n) vra | n G ^} U {(a2 + n)7!-a \ n E Z}. The operator Xa^^aj has the 
complete set of eigenfunctions 

fi"i\x) = ^J— / cos{ak) exp[— i (ai + n) tt sin(aA;) + ikx]dk (6.13) 



2 V TT 

1 Ici 



2a 



f!i^2\^) = 2\l \J ~^ J j \/l cos(aA;)| exp[— i (0:2 + n) vr sin(aA;) + ikx]dk (6.14) 

\ a 2a / 

in £2(0^)- For the umbral map we can draw the following conclusions, 
operator: y ^ai,a2 

spectrum: H {{ai + n)'Ka\nEZ, £ = 1,2} 

eigenfunctions: fx{y) = S{y - X) fi,j\ fl^^^ . 

Of course, in this case we have no method to calculate eigenfunctions of Xoi,q2 directly 
from the generalized eigenfunctions 5{y — A) of the Schrodinger operator y. 

Slightly more complicated is the case of the operator x^. Following our treatment of 
the operator x itself, a set of two selfadjoint extensions of the operator —d'^/ dp^ determines 
a selfadjoint extension of x^. The domains of selfadjoint extensions of —d^/dp^ in L^^^ 
have the form 



dx , , d^x 

— abs. cont., —7 
dp dp^ 



T^b.c. = {x G L]^!^ I X differentiable, — — abs. cont., -—r G L^^u, h.c.} (6.15) 



where h.c. stands for a certain choice of boundary conditions, like x(^l/'^) = = x{^/(^) 
(see 0] for other choices). 

Example. Let us consider the equation aapQ = ki/jq where a. = d/dy + y is the annihilation 
operator for the one-dimensional harmonic oscillator, and k G (D. The umbral map 
replaces a by Q + x. Choosing for x a selfadjoint extension, we have to consider 

(Q + Xai^aJ V'O = K V'o . (6.16) 
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We write the Fourier transform of ipQ as ^E'ol^) = y \ cos{ak)\ Xoip) with p given by 
( |675|) , separately on [— 7r/(2a), 7r/(2a)] and [— vr/a, — 7r/(2a)] U [7r/(2a), vr/a]. Now (|6.16|) 
translates on both subsets of [—vr/a, vr/a] to {p + d/dp)xo = /^Xo with the solution 
Xo{p) = C exp{Kp — p^ /2) where C is a constant. For C 7^ one finds Xo ^ T^a with 
a = — i/t/(vra). This restricts k since a G [0, 1). Furthermore, ai = 02 = Application 
of the (umbral-mapped) creation operator —Q + ^a,a to t/'q leaves V^^^ since pxo(p) is not 
in Va. The algebraic construction of the eigenfunctions for the harmonic oscillator there- 
fore does not survive after the umbral mapping. The problem actually appears already 
in rewriting the Hamiltonian as 

Ha ■=\{-Q^ = \{-Q + ^a,a){Q + ^a,a) + \[Q,^a,a\ ■ (6.17) 

The point is that [(5,Xa,o] = I does not hold on the domain of Xq,,^. As a consequence, 
there is no simple relation between the spectra of Ha and | {—Q + ia,a){Q + ^a,a)- 

In the way described above, the eigenvalue problem for a selfadjoint extension of the 
Hamiltonian H = {—Q^ + x^)/2 reduces in p-space to (twice) the eigenvalue problem 
for the Hamiltonian of the ordinary harmonic oscillator restricted to the finite interval 
[— 1/a, 1/a] with the respective boundary conditions. A choice among the many different 
selfadjoint extensions of H should be determined by the specification of the physical 
system (on the lattice) which we intend to describe. It is not obvious for us, however, 
what a natural choice could be. D 

An interesting aspect of the representation of the canonical commutation relations 
considered in this section is the fact that it is solely composed of the two operators x and 
Q which both receive a physical meaning if we interprete i2{aZ) as the space of functions 
on a (physical) space lattice, x is the position operator and —iQ the natural candidate 



for the momentum operator (see also [T^). This is the basis for a discrete version of 
quantum mechanics. Whereas ordinary quantum mechanics has a continuous position 
space, discrete quantum mechanics lives on a lattice. Quantum mechanical models on a 
lattice should then be modelled with the selfadjoint operators x and — iQ. These satisfy 
commutation relations which are different from the canonical ones. Still missing is, how- 
ever, a general recipe to quantize a (discrete) mechanical system, analogous to canonical 
quantization. But what is the meaning of the representation given by x and —iQ ? Ba- 
sically it just offers us a way to get, apparently, close to the results of ordinary quantum 
mechanics within the framework of discrete quantum mechanics. That this representation 
is not equivalent to the Schrodinger representation means that, within the framework of 
discrete quantum mechanics, we cannot reproduce ordinary quantum mechanics rigor- 
ously, at least not in the way attempted in this section. In fact, we have found rather 
drastic deviations, in particular a kind of spectrum doubling, a familiar problem in lattice 
field theories & ITDI . 
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7 Conclusions 



In this paper we have pointed out that there is an apparently widely unknown math- 
ematical scheme, called umbral calculus, behind recent work 0, |^ on discretization of 
differential equations and physical continuum models. Using several examples we have 
discussed its prospects and shortcomings. By choosing delta operators different from 
those used in these papers, alternative discretizations can be obtained. They have not 
been worked out in detail yet. 

Discretization of a continuum theory breaks the continuous space-time symmetries 
which play a crucial role in (non-gravitational) quantum field theory. There have been 
attempts to find a discrete analogue of space-time symmetries for lattice theories such 
that essential features of the continuum group structures are maintained. Discretizations 
of Lorentz transformations were considered in fl^, for example. In [|l^ the Poincare 



group acts on an ensemble of lattices (see also [|I^ for a related point of view). Umbral 
calculus offers a different way to implement symmetries on lattices (see also P]). 

Umbral calculus provides us with certain classes of representations of the canonical 
commutation relations. It is therefore of potential interest for quantum mechanics and 
quantum field theory. Among the variety of umbral maps which we have at our disposal, 
the one determined by ( ^.4|) with the central difference operator is of special interest (see 
also [^]). In this case we have a representation of the canonical commutation relations 
constructed from the position and the momentum operator on a lattice. This suggested 
a kind of embedding of ordinary quantum mechanics into a formalism for quantum me- 
chanics on a lattice and thus a discretization of quantum mechanical systems which is 



different from conventional ones (see |19|, for example). The representation of the canoni- 
cal commutation relations obtained in this way is, however, not unitarily equivalent to the 
Schrodinger representation. As a consequence, the image of ordinary quantum mechanics 
under the umbral map cannot reproduce the results of the former rigorously. We revealed 
a kind of spectrum doubling similar to what is known in lattice field theories. This may be 
regarded as a negative feature. In any case, we believe that this is an interesting example 
of a representation of the canonical commutation relations by selfadjoint operators which 
is not equivalent to the Schrodinger representation. Furthermore, our analysis sheds some 
light on the work in |^ where this representation has been used. The umbral framework 
yields many more examples, of course, which can be analyzed analogously to the example 
which we selected in section 

Acknowledgment. F. M.-H. is grateful to Klaus Baumann for some helpful discussions. 



A On the *-product associated with the forward dis- 
crete derivative delta operator 

Let X be the canonical coordinate function on a lattice with spacings a. A function f{x) 
for which the finite difference analogue of the Taylor series expansion (Gregory-Newton 
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formula) exists can be written as follows, 

oo oo 

f(x) = J2 Ux- = J2F,x^'^ (A.l) 

n=0 k=0 

where we have used 

n 

= ^ ^(n,fc)a"-^x('=) . (A.2) 

A:=0 

The coefficients S{n, k) are the Stirling numbers of the second kind {S{n, 0) = when 
n > 0, S{n,n) = 1). The coefficients Fk in ( [A.l| ) are given by 



Fk=T. fk+mS{k + m,k)a^ . (A.3) 



m=0 



The equation ( |A.1| ) can also be expressed as 

/(x) = F{x) (A.4) 

where 

oo 

F(x) = 5:FfcX^ (A.5) 

fc=0 

The *-product of two functions of x is then given by 

oo 

(/ * h){x) = {FH){x) = FkHex^'-"'^ (A.6) 

k,e=o 

and the rhs can be written as a power series in x with the help of 

n 

x^'^) = ^ s(n,A;)a"-'=x'= (A.7) 

k=0 

where s{n, k) are the Stirling numbers of the first kind. 
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